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Abstract. 

' This is a brief review of applications of extended gravity theories to cosmology, in particular to 

the dark energy problem. Genetically extensions of gravity action involve higher derivative terms, 
■ which can result in ghosts and instabilities. We consider three ways to circumvent this: Chern- 

Simons terms, first order variational principle and nonlocality. We consider some recent cosmolog- 
ical applications of these three classes of modified gravity models. The viable parameter spaces can 
be very efficiently bounded by taking into account cosmological constraints from all epochs in ad- 
dition to bounds from Solar system tests and stability considerations. We make some new remarks 
concerning so called algebraic scalar-tensor theories, biscalar reformulation of nonlocal actions in- 
volving the inverse d'Alembertian, and a possible covariant formulation holographic cosmology 
with nonperturbative gravity. 
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1. INTRODUCTION 



General theory of relativity (GR) has to be extended at high energy scales to resolve its 
singularities, the hierarchy problem, quantization and unification with other theories of 
fundamental interactions. The dark energy energy problem [1] has raised speculations 
about the possibility of an infrared modification of GR too. Alternative theories of 
gravitation have been around since Nordstroms scalar gravitation [2], but no evidence 
for departure from the classical predictions of GR have been observed. One may then 
contemplate if the acceleration of the universe (and perhaps also the need to invoke dark 
matter) could be a sign of a need to modify GR. 

In Uzan's review of possible physics behind the cosmic acceleration [3], a classifica- 
tion of different approaches was based on which kind of new fields they introduce and 
how these new fields couple to the matter and the metric. Generically, modifying gravity 
is equivalent to invoking a new matter field. In the simplest case of a Brans-Dicke -like 
theory [4] it is a scalar field but other types of fields can appear as well, for example 
tensors in bimetric theories. Usually also several fields come about, like in low-energy 
limit of string theory featuring typically form fields, the dilaton scalar field and other 
moduli. 

In the present discussion we will focus on models which avoid introducing new fields 
while modifying gravitational interactions. The motivation for this is that, as will be 
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explained in the next section, rather generically the new fields would actually behave 
pathologically and thus actually render the starting point invalid. It is interesting to 
contemplate how well such specific classes of gravity theories are constrained, since 
in principle all other theories can be understood as new matter in disguise (as far as 
only gravity is concerned). There is also some appeal in the attempts to deal away dark 
energy without invoking any fields, though of course no fundamental principle excludes 
the possibility of the dark sector being arbitrarily complicated: indeed it could involve as 
many fields as the visible sector described by the Standard model - or even vast number 
of more fields [5]. 

This starting point complements the reviews of modified gravity cosmology [6, 7, 
8, 9, 10]. The plan of the presentation is as follows. In Section 2 we present simple 
field theory considerations which provide the basis for our discussion. The lesson is 
that unless new fields were introduced, only specific classes of effective field theories 
are viable. We then look at this within the scalar-tensor framework. Then three classes 
of models are described and a brief summary of their cosmological applications to the 
dark energy problem is provided: the Gauss-Bonnet model in Section 3, the first order 
formulation of modified gravity in Section 4 and nonlocal gravity models in Section 5. 
We conclude in Section 6. 

2. EXTENSIONS THE GRAVITY ACTION 

Consider a generic Lagrangian involving derivatives to order N, L(q,q, . . .,qW). The 
canonical phase space then has 2N coordinates 1 , and they may be defined as 




One might use the equation of motion to solve = A(Qi, . . . ,Qn,Pn) in terms of 
the other variables. However, the other degrees of freedom remain unconstrained and 
contribute to the Hamiltonian with indefinite sign: 

N 

H = -L = P1Q2 + P2Q3 + • • ■ + Pn-iQn+PnA-L(Q u QnA) • (2) 

(=1 

Thus the Hamiltonian is unbounded from below. This result is generic, unless additional 
constrains are introduced [11], higher derivative theories lack a minimum energy state. 
This may be called the Ostrogradski theorem. See Woodard's lecture [12] for a lucid 
presentation of the argument and discussion within the framework of generalized gravity 
theories. 

Below we will consider a consider prototype of generalized gravity, scalar-tensor 
theory, from the point of view of eliminating or constraining the additional degree of 
freedom. There is a one-parameter family of classes of scalar-tensor theories in which 



We assume the Lagrangian depends nondegenerately upon . 
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the scalar field is given by an algebraic function of the curvature and the trace of matter 
content. The family of these theories provides a new perspective to the nonlinear gravity 
theories, since it includes both their metric and Palatini formulations as special cases. It 
is then possible to continuously interpolate between these formulations and avoid some 
pathologies. To end this section, we resume in 2.2 the possible approaches to modify 
gravity avoiding the Ostrogradskian catastrophy. 



2.1. A scalar- tensor viewpoint 

We consider the scalar-tensor action in n dimensions 



1 1 . ., 2 



(3) 



where ^f( m ) is the matter Lagrangian and the corresponding stress energy tensor is 
defined as 



The field equation for the metric follows as 



(V M 0)(V v 0)--(d0) 2 g M v 



F(0)G MV = (V„V v -g MV [D)F(0) + a>(0) 

- V($)g^ v + W($)T$. 
The scalar field obeys the equation of motion 

+ l -CD'md<l>f - + l -F'{^)R = -W(0)J*f ( „ 
From these equations one may infer that the matter sector obeys 

VTW = (VlogW(*)) (^ W -#) , 
expressing the co variant energy-momentum conservation [13]. 



(4) 



(5) 



(6) 



(7) 



2.7.7. The algebraic family 

By combining the trace of the field equation (5) and the generalized Klein-Gordon 
equation for the scalar field (6), the scalar curvature can be expressed as 



a 



(1-\)F-Un-1) F - 
2 2 CO 



R 



(l-^)co+ l -(n-l)co'^-(n-l)F" 
2 2 CO 



nV-{n-l)^-V + WT {m) + {n-l) ^ WSf {m fi) 

CO y ' CO 
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We ease notation, we will occasionally drop the explicit -dependence from the func- 
tions as we did here. If the first square bracket term in the LHS vanishes, there is an 
algebraic relation between the curvature, the matter variables and the scalar field. Oth- 
erwise such a unique relation does not exist, though it is satisfied by any theory in the 
dynamical situations where the derivatives of the field happen to vanish. If the relation 
holds universally within the theory, we call it an algebraic scalar-tensor theory. 

To analyse their properties, it is convenient to set F(<j>) = 0. We can do this without 
any essential loss of generality, since as long as F((j)) is invertible, the change is equiv- 
alent to a redefinition of 0), V and W. Now it becomes easy to show that the algebraic 
theories allow only a certain form of the kinetic term. In particular, the requirement is 
that 

where Q.^ is an arbitrary constant (stemming from an integration constant), and the 
constant Q. n depends on the dimensionality as 

n-l 

a n = — -. (io) 

n — 2 

The constant Q.^ G R is allowed to take arbitrary values. Thus the algebraic theories 
constitute a one-parameter family of classes of theories in the space of all possible scalar- 
tensor theories. 

The nonlinear gravity models, the so called f(R) gravities, are specific examples 
of these algebraic theories. In the metric formulation of nonlinear gravity, the kinetic 
term is absent and thus this corresponds to the limit of infinite Q.a- The field is still 
dynamical, because it is nonminimally coupled to gravity. However, it can be expressed 
a function of a R alone, § = §(R). It has been argued that indeed in this particular case 
one may succeed in avoiding the Ostrograskian instabilities while increasing the order of 
derivatives [12]. However, in the present discussion we do not go into details of this case. 
Extensive studies are found in the literature (for a list of reviews, see the introduction). 

The case corresponding to the limit of vanishing Q.a is nothing but the "co = — 3/2 
Brans-Dicke theory" and it is in many regards the opposite of the metric f(R) models. 
This is equivalent 2 to the first order formulation of nonlinear gravity (i.e. Palatini-/ (R) 
models) which we discuss in more detail in section 4. In such a scalar-tensor theory then 
the field is not dynamical, since the kinetic term exactly cancels the effect of nonminimal 
gravity coupling. One may express the field as a function of T alone, = (j)(T). 

As one interpolates between the extremes £Ia = and Q.^ — > °°, the scalar field always 
depends on the both the trace of matter and the Ricci curvature. However, these theories 
become equivalent to metric nonlinear gravity in vacuum or in the presence of conformal 
matter; and they are equivalent to Palatini nonlinear gravity in Ricci-flat spacetimes. 
The algebraic scalar-tensor theories are conformally coupled scalar models, whose the 
gravitational constant is given by Q.^. 



2 Or slightly more general, since the potential can be arbitrary and does not have to be obtained from a 
particular function f(R). 
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2.2. Three approaches 



There are three frameworks of extending gravity theories without introducing new 
degrees of freedom, in particular avoiding ghosts. 

• Gauss-Bonnet models. The Gauss-Bonnet invariant G = R 2 — 4R^ V R^ V + 
Rp.vpoR^ vpa combines the quadratic invariants in such a way, unique in four 
dimensions, that the higher derivative contributions cancel. Models based on this 
term could thus avoid higher derivatives. 

• Palatini models. In the Palatini framework one regards connection and the metric as 
independent variables. This means working on first order formalism instead of the 
second order one, or in Hamiltonian instead of Lagrangian formulation. Then one 
can consider extended action without necessarily increasing the order of the theory. 

• Nonlocal models. For an action involving infinite number of derivatives, the 
Ostrogradski-type argument above can break down as a truncation to cannot be 
performed at any finite order to construct the canonical momentum. Such nonlocal 
models can modify the gravitational interaction without new propagating modes. 

Gravitational alternatives to dark energy belonging to each of these frameworks have re- 
ceived considerable attentation in the literature during few recent years. In the following 
we review concisely what has been learned about their cosmological implications, not 
aiming at an exhaustive account, but a summary of the main features of some presenta- 
tive examples. 



3. GAUSS-BONNET MODELS 

The action for the system to consider is 

v 

(ID 




M 2 v 

- . jR _Z(v^,) 2 -y(0)-/(^,)G+L„ 



2 2 

where 7 is a constant and is the coupling which appears in the one-loop corrected 
string effective Einstein frame action [14, 15]. The Gauss-Bonnet term G is topological 
and thus doesn't alone contribute to the field equations. Therefore we cheat by including 
the field 0, which can be identified with the string theory dilaton or a modulus. We 
parameterize the coupling and the potential as 

V(<j)) = Vo e x ^ M \ f(<j>) = /^/(v^m) (12) 

The nonperturbative effects from gaugino condensation or instanton can result an expo- 
nential potential. An exponential field dependence of the coupling is known to occur in 
supergravity actions. These forms are thus theoretically well motivated. We will see that 
the exponential case happens to have some beneficial features for phenomenology too. 
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TABLE 1. The critical points in the system. 



Fixed point 


X 


y 


J" 


w eff 


A 











Wm 


B 


±3 








1 


C 





3 





-1 


D 


A/2 


3-A 2 /4 





A 2 /6-l 


E 


3(l+w m )/A 


9(1-m4)/A 2 





w m 


F 


3(1 +w m )/a 







w m 


G 


^0 





^0 


w 



3.1. Gauss-Bonnet models: Cosmology 

Though the severe constraints on Gauss-Bonnet coupling of dark energy were men- 
tioned in [16, 17], the possibility of Gauss-Bonnet interaction of dark energy has been 
considered in several studies [18, 19, 20, 21]. Here, in particular we discuss the scenario 
proposed in [22], see also [23]. Other models include f(G) gravity [24], generalized 
couplings [25], vector couplings [26] and brane- world models [27]. 

We consider the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

ds 2 = -dt 2 + a 2 (t)dx 2 (13) 

and refer to H = d/a as the Hubble parameter, an overdot meaning derivative with 
respect to t. The equations of motion for the scale factor and the scalar field are then 

3M 2 H 2 = ^<j> 2 + a 2 V(<j))+24a- 2 H 3 f / (<t))<t) + p m , (14) 

7(0 + 2H<j)) +a 2 V'(§)+2Aa- 2 f\<l>)H 2 H = 0. (15) 

For the purposes of phase space analysis, it is convenient to introduce the dimensionless 
variables defined as 

O " K2pm r~ K * v-K- 2V ^ 



3H 2 ' y/2H' 7 H 2 



H = 8^0^/(0), Weff = -11L-1. (16) 

The equations (14,15) may then be rewritten in terms of such variables [21], and the 
fixed points of the dynamical system that we find corresponding to the exponential 
model (12) are then described in the Table 1. The stability properties of these fixed 
points are summarized in Table 2, for the details see Ref.[28]. The dark energy scenario 
is basically the transition from the scaling fixed point E to the accelerated fixed point G, 
as we describe below. 

The exponential quintessence is known to mimic the background dominating com- 
ponent equation of state. This so called scaling solution is an attractor, and thus the 
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TABLE 2. Description of the critical points in the system and a 
comment on their stability or existence. 



Fixed point 


stable? 


description 


A 


no 


matter domination 


B 


no 


kination 


C 


if a > X 


de Sitter 


D 


doesn't exist for X > 2\/3 


scalar dominated 


E 


saddle point if a > X scaling 


F 


no 


G-B scaling 


G 


no 


G-B domination 



3.5 




FIGURE 1. Left panel depicts the fractional energy densities for matter, i2,„ (dash-dotted line), the 
scalar field Q.^ (dashed line) and Gauss-Bonnet term, Q.f (dotted line). The solid line is the total equation 
of state w e ff. and the transient phantom era is caused by the dynamics of the coupling. Right panel is a 
portrait of the phase space. Here X = 4 and a = 20. 



evolution of the field is widely independent of the initial conditions [29]. To accelerate 
the universe however, the field should be tossed of the attractor. Now this is naturally 
done by the Gauss-Bonnet coupling, if a > X. Given this simple condition, the scaling 
attractor becomes a saddle point from which the universe enters into the asymptotically 
de Sitter space. Possible transient effects include phantom-like expansion with w < — 1. 
The full evolution of the universe is shown in the figure 1 . 

It is pleasing that not only the form of the model but also the numerical values relevant 
for this scenario can be motivated by low energy string actions, since the parameters 
we need for the present acceleration in (12) have natural scales. The dimensionless 
parameters a and X are of the order one, and due to the shift symmetry of the potential, 
one can always consider Vq ~ M 4 and fo ~ M~ 4 . The initial conditions of the field 
are (almost) arbitrary. Similar exploitation of the exponential form can be deviced by 
adjusting the form of the potential of the field or disforming its coupling to the metric 
(when in fact a similar transition mechanism emerges [30]). 

The Solar system limit of these models is problematical, as already observed by 
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Esposito-Farese [16, 17]. The Gauss-Bonnet term behaves as r 4 , and it is then difficult 
to neglect its small scale effects. Recently ways out have been considered, by recon- 
structing couplings that reproduce GR at small scales [31] or applying the chameleon 
mechanism [32]. At the cosmological scales new effects may occur. The evolution equa- 
tion for the matter overdensity 8 is given by 

S + 2H8 = 47cG,p8. (17) 

The effective gravitational constant G* seen by the matter inhomogeneities depends 
rather non-trivially on the evolution of the background quantities, and in terms of our 
dimensionless variables defined in Eq.(16). It can then be written as 

r -* 4 + M 2 (l+£) 2 + * 2 [2(l+£)(M-l)+y] (m 
* * 2 [4 + J u(5 J u-8)]-Ai 2 [6(l+£)( i u-l)+y] ^ ; 

In fact, it is possible that the expression for the denominator vanishes, and the effective 
gravitational constant diverges. This is linked with an appearance of a ghost [28]. 
Such peculiar instabilities exist both for tensor and scalar perturbations in one-loop 
string cosmology. Some of these instabilities were first discovered in the pre-big bang 
cosmology [33]. Recent considerations of stability include [34, 35, 36, 37]. 



4. PALATINI MODELS 

Consider gravity theories represented by the action 3 



S = (Fx 



(19) 



Here 0,... are some matter fields. In the Palatini approach one lets the torsionless 
connection n? vary independently of the metric. The Ricci tensor is then constructed 
solely from this connection, 



R 



fa _ fa 
1 nv,a 1 )ia,v 



+ 1 rvl 1 



aA A Atv 1 /jA 1 av ■ 



(20) 



In the following we will call R = g^ v R^ v . The field equations which follow from 
extremization of the action, Eq.(19), with respect to metric variations, can be written 

as 

FB$-\f8$ = T$, (21) 



3 More general functions f{g llv ,R^v) can be considered too [38, 39, 40], but we restrict here to the sim- 
plest case. In addition, if matter was nonminimally coupled to the independent connection we would 
be lead to more general models characterized by nonzero hypermomentum [41]. For our purposes of 
demonstrating cosmological phenomenology the simplest case is enough, however to address the criti- 
cisms concerning the validity of the averaging procedure [42], existence of simplest star solutions [43] 
and well-posedness of the Cauchy problem [44], one may be forced into a more general framework. 
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where we have defined F = df/dR. In GR, / = K 2 (R - 2A), so F = K 2 . By varying the 
action with respect to rjj^, one gets the condition 



-88 



0, 



implying that this connection is compatible with the conformal metric 



(22) 



(23) 



This connection governs how the tensor v appearing in the action settles itself in order 
to minimize the action, but it turns out that the metric compatible connection determines 
the geodesies that freely falling particles follow, since the energy momentum defined 
as (4) is conserved according to this connection [13], V^T^ 1 = 0, whereas in general 

VnTy 7^ 0. Therefore one may interpret the Levi-Civita connection as the gravitational 
field as in Ref.[45]. The trace of the field equations allows us to solve R as an algebraic 
function of the matter trace T = g^T^. This central relation is 



FR — 2f=T. 

Written in the form of GR plus correction terms, the field equations read: 



(24) 



<%(g) = Tt + (l-F)R>i( g )-—(WF)(V v F) 



(f-R) + (\-j)nF + ^(dF) 2 



8». 



(25) 



Since the corrections can be expressed as functions of the matter trace, one can view 
Eq.(25) as GR with nonstandard matter couplings 4 . The whole RHS may be regarded 
as an effective matter energy-momentum tensor. In vacuum it reduces to a cosmological 
constant [47]. This is also the case in the presence of conformal matter, i.e. if T = 0. 



4.1. Palatini models: Cosmology 

In Palatini-f(R) cosmology [48, 49], the Friedmann equation obtained from (25) can 
be written as 

1F^ 2 



6F(H + --) =p+3p-f. (26) 



If w is constant this may be rewritten as [50] 



6H 2 = 1 (1+3mQ*F-3(1+mO/ 

\ l 2^ ' F{RF" —F) J 



4 An unusual approach with possibly similar consequences is to consider f(R, T) Lagrangians [46]. 
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We have expressed the Hubble rate as a function of R, which we in turn may solve 
from the trace equation (24). Thus it is straighforward to find the scale factor evolution 
given the matter content. It turns then out that these models can generate a viable 
sequence 5 of radiation dominated, matter dominated and accelerating era matching with 
the constraints. This has been shown for various parameterizations of the function f(R), 
most often with some power-law forms (with one, two or three powers of R), but also 
with square-root, logarithmic and exponential forms for the curvature correction terms 
[53, 50, 54, 55, 56]. As an example we show, in the left panel of Fig. 2, the constraints 
arising from fitting the the last scattering distance data to the prediction of the model 
defined by a monomial correction to the Einstein-Hilbert action. 

At the Solar system scales, it seems that these models can correctly reproduce the GR 
phenomenology. Since the gravity correction in vacuum reduce to GR with a cosmolog- 
ical constant, and the magnitude of the effective cosmological constant is given by the 
scale of acceleration today, the models can pass the local gravity tests. 

The problems occur at the level of inhomogeneous perturbations. For the background, 
the additional derivative terms in the field equations due to nonlinearity in f(R) can play 
the role of an effective smooth dark energy, but the spatial gradients of these extra terms 
cannot be neglected for cosmological perturbations [57, 58] which then assume unusual 
behaviour that is at odds with the observed distribution of galaxies [59] and with the 
cosmic microwave background spectrum [60], even if the nonlinear part of the action 
is exponentially suppressed at late times [61]. In quantitative terms, parameterizing the 
f(R) — R ~ RP, data analysis constrains |j8| < 10~ 5 [57, 60]. We show the confidence 
contours arising from matching with SDSS data [62] in the right panel of Fig. 2. 

However, there is way out of the problem. If we consider a dark matter component 
which is not entirely cold, but may have finite sonic, anisotropic or entropic stresses at 
large scales the effects of the modified gravity on the matter power spectrum may be 
partially or completely eliminated. It is then possible to obtain exactly the same linear 
matter growth as in a universe with smooth dark energy producing the same background 
expansion. Such Palatini models with generalized dark energy may be called f(R)GDM 
models instead of f(R)CDM models which do not are not viable alternatives to the 
ACDM models. The evolution of the dark matter overdensity is, in general coupled to 
the evolution of the anisotropic stress IT as 



The exact form of the five dimensionless coefficients appearing in this equation can 
be read from the appendix of the Ref.[63], where it was shown that any f(R)GDM 
satisfying the relation (the sound speed c 2 can vanish as for CDM) 



A simple quadratic term can also generate a bouncing cosmology [51], showing these models might 
have relevance as an effective description of loop quantum cosmology [52]. 



'8 = DiH8 + (D 2 H 2 + k 2 D k ) 8+P 1 HTl + P 2 H 2 Y\. 



(28) 




(29) 
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FIGURE 2. Constraints on the power-law Palatini model f(R) = R + aRp. Left panel shows the back- 
ground constraints arising from from fitting the CMB shift parameter. Right panel shows the perturbation 
constraints ensuing from fitting to the observed matter power spectrum. The contours correspond to 68, 
90 and 99 % confidence levels. We allowed generous variation of the spectral index ns = 1 ± 0.2. Note 
that the j3-axis of the right panel is scaled by 10~ 5 . 



escapes the matter power spectrum constraints completely. However, due to presence 
of anisotropic stress now both in the matter and in the gravity sectors [64, 65], the 
gravitational potentials will respond differently to the clumps of matter. These potentials 
affect the lensing experiments and also the CMB via the integrated Sachs- Wolfe effect 
[66]. Thus one may distisguish f(R)GDM from ACDM by exploiting these data. The 
former models have not yet been confronted with the data. A recent data analysis 
using a general parameterization [67] favors standard growth rate but nonzero effective 
anisotropic stress, which is the generic /(7?)GDM prediction. 



Thus the model is specified by the infinite number of coefficients c n . One could also 

write them as c n = c n /M^ n+1 \ where c n are dimensionless and if roughly of order one, 
M* is the mass scale at which the higher derivative terms in the action become impor- 
tant. To avoid the generic higher-derivative problems, the coefficients c n must remain 
non-vanishing to arbitrarily high orders. This string-inspired form of nonperturbative 
gravity has been also shown to be asymptotically free and thereby potentially describe 
bouncing cosmologies [68, 69]. Similar models may be also motivated by quantum loop 



5. NONLOCAL MODELS 



A class of nonlocal gravity models can be defined by the action 




(30) 
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corrections to gravity. We consider then the following simple example 6 of a nonlocal 
action [74, 75] 

S = J d A x^[R{\+f(n- l R))+2K 2 J? m ], (31) 

where / can be called the nonlocal distortion function. If we consider analytic functions 

/(*). " 

f(x) = J £f i x\ (32) 
(=0 

then the actions can be shown to reduce, at least patchwise, within the class (30). Con- 
sider the simple case f(x) = fix. Such correction was proposed to bound the Euclidean 
action and its cosmological effects were considered already in Ref.[76]. This model is 
equivalent to the series specified by 

* = — *E(^*)T <33) 

One should note here that the scale must be such that —DR <C M 2 R. In a similar way 
it is straightforward to see that with the same caveat, any analytic / can be presented via 
such a power series (the coefficients c ; are then given by a double sum whose explicit 
form is not needed here). 

Introducing the Lagrange multiplier t, to rename the inverse d'Alembertian as 0, we 
may rewrite the action in the form 

S = Jd 4 x^[R(\+f($)) + t;(n$-R)+2K 2 J? m \. (34) 

This is a peculiar kind of double- scalar-tensor theory where massless, nonminimally 
coupled scalars have a kinetic-type interaction resembling of the nonlinear sigma model. 
If we then define 

= (35) 

integrate partially and neglect the boundary terms as usual, it it becomes transparent that 
we may consider a biscalar-tensor theory, with only one of the scalar degrees of freedom 
is nonminimally coupled to gravity as 



S = |j 4 x v ^[(l+^-/(0)((90) 2 + V a W O! 0+2jc 2 ^ w ] . 



(36) 



However, like in reducing the action into the form (30), we should note some caveats in 
reducing it into the form (36). Namely, when getting rid of the inverse d'Alembertian in 
step (34), we are introducing a second order constraint into the Lagrangian. In general, 
this is not legitimate, as pointed out by Koshelev [77]. The solution to the constraint 
equation D0 = R is unique only up to the most general harmonic function <p for which 



6 For cosmology of nonlocal scalar field models motivated by string field theory see for example [70, 71, 
72,73]. 
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□ <p = in a given background. Thus the local form of the model seems to introduce 
a new degree of freedom which the original theory lacks: the Lagrangian multiplier 
becomes a dynamical variable. From the practical point of view of finding classical so- 
lutions the relevant question is whether one can recover the solutions of the original 
theory by just choosing the initial conditions suitably, i.e. choosing the correct function 
(p = 0. As a particular case, consider the limit / — > 0. The original theory obviously re- 
duces to GR, while the spurious degrees of freedom remain in the action (36): however, 
it is clear that the latter are not actually there if we impose the correct boundary con- 
dition for the equation of motion 7 . This is the only case where the non-equivalence has 
been explicitly shown [77], but a more general treatment of the issue would be worth- 
while, especially since the local form (36) (with the prescription for initial conditions 
understood) is considerably more tractable than the manifestly nonlocal starting point 
(31). A rigorous mathematical analysis of the possible equivalence 8 is outside the scope 
of present discussion, and thus one should regard some the following results strictly to 
apply only to the local variant (36) of the model (31). 

5.1. Nonlocal models: Cosmology 

The reasons to suspect a correction involving R/D as a culprit for dark energy can 
be easily seen. This term is dimensionless. Thus, one may construct the distortion 
function without introducing any dimensional functions. It appears as a prefactor for 
the Newton's constant, and explain weakening of gravity at cosmological scales, the 
distortion function should be about minus few dozens of percents today. Thus no new 
scale needs to be introduced, in particular the hierarchy between the observed magnitude 
of dark energy density and the Planck scale does not to have to prescribed into the 
action. Another side of the problem is then why the correction is important at the present 
stage, and in this model the explanation is that during the radiation dominated era the 
Ricci scalar is vanishing. So the evolution of the distortion function and the ensuing 
acceleration is a consequence of the universe entering into the matter dominated stage. 
-The possible problems are also easily seen. Namely, it is not clear that the distortion 
effects would be confined solely to the cosmological scales and the question arises why 
there aren't observed elsewhere. 

Analytic and numerical solutions for the background expansion in this model in a 
realistic universe were found in [83], see also [84, 85, 86]. A method to reconstruct 
the model generating an arbitrary background history was also constructed for both the 
scalar-tensor form and the original model [87]. Another way to proceed is to specify 
a function / and compute its predictions for the FRW metric. For simplicity, consider 
the monomial f(x) = fnx" (no summation). We may then schematically summarize the 



In this simplest (and trivial) case, the equation of motion for is UN* = and the correct boundary 
condition reproducing the original model is *P = 0. 

8 Recent localization considerations [78, 79, 80] have for example employed reformulation of the nonlocal 
system as diffusion problem [81, 82]. 
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possible evolution histories as follows: 



n 



n>0 



Nonlocal effect 



n<0 



(-i)"/»>o 




Matter domination 



Slows down expansion 



Acceleration 



Singularity 



Thus, with positive the power n and the sign /„ equal to 9 — ( — 1)'\ one finds dark energy 
solutions. This scenario leads to a sudden future singularity, which might be avoided be 
regularization of the operator or reshaping the function /. Though this was not in the 
focus of the present study, n < may lead to a decaying effect which, with large enough 
scale /„, can influence the evolution of the early universe. 

The linear inhomogeneity field of matter perturbed around the metric (13) evolves as 



The brackets in the last expression enclose the effective gravity coupling. We see it 
is time-dependent but scale-invariant. We then do get nontrivial modifications of the 
growth rate, but don't find the problematic effects of a nonzero sound speed like in 
the first order gravity models of section 4.1. The time-dependence of the gravitational 
constant can in principle be used to distinguish a nonlocal origin of acceleration, when 
the background expansion is identical to some local model. 

The perturbations about the Schwarzchild geometry for the model (32) can be shown 
to imply the post-Newtonian parameters 



Therefore it is possible to constrain the constant and the linear parts of the coupling 
from Solar System experiments. However, the constant fo may always be absorbed into 
a redefinition of the jc — > ?c/(l +/o), and so one may set fo = 0. We then get the strict 
bound —5.8 ■ 10 6 < fi < 5.7 • 10 6 . Deviations from GR ensuing from the higher order 
terms in the expansion for f(<j>) are of the second order in perturbations, ~ (G^M/r) 2 , 
and thus seem undetectable within present experimental accuracy. However, we have 
only shown that such solutions exist. Whether the assumption that we may expand about 
the Schwarzchild geometry is realistic should be assessed with more care. Properties of 
exact spherically symmetric solutions were studied in [89] . 



5.1.1. Nonperturbative gravity and holographic cosmology 

It is interesting to remark that nonperturbative string-inspired (30) or quantum loop 
corrected (31) gravity models could provide a co variant realization of the idea of holo- 



y In fact this sign seems necessary to avoid a ghost. 
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[88] 




(37) 




(38) 



graphic cosmology. The latter is indeed some manifestation of nonlocal gravity. As we 
explained above, the nonlocal distortion function can automatically feature effects at 
the right energy scale, thus alleviating the cosmological fine-tuning problems. A similar 
promising feature can be found in the holographic dark energy models. The holographic 
principle, familiar from black holes [90] and quantum gravity [91], is postulated to set 
an upper bound on the entropy of the universe [92] . The quantum zero point energy of 
the region of size L should not then exceed the mass of a black hole of the same size, 
L 3 Pa < L/k 4 . This then implies dark energy density scaling as L 2 , and if L is propor- 
tional to the Hubble horizon L ~ l/H, the fact that dark energy density today happens 
to be of the observed magnitude becomes almost tautological. However, the details do 
not work correctly unless the horizon is given by the future event horizon [93], 



This is the boundary of the volume an observer at time t may eventually see. At this 
level, we only have an arguably strange modification of the Friedmann equation which 
lacks a covariant justification. However, it is not unconceivable that similar cosmologies 
may be derivable from a nonlocal theory of the form (30) or its variants. In particular, 
if one takes the action (31) and considers the inverse d'Alembertian to correspond to, 
instead of the retarded Green function, the advanced one, then future integrals like those 
in the expression for the event horizon (39), appear into the effective field equations. For 
example, in FLRW background we have 



Producing exactly the modification (39) seems difficult 10 but is probably unnecessary, 
as the main aim would be to produce dynamical dark energy whose magnitude would 
depend upon some suitably defined horizon size. 



Modified gravity often contains pathological degrees of freedom originating from higher 
derivatives or extra fields. Avoiding introduction of new propagating modes is possible 
in some special cases, whose cosmological implications were considered here. In the 
framework of scalar-tensor theories, a special one-parameter family of theories is formed 
by those in which the scalar field is an algebraic function of the Ricci scalar and the trace 
of the matter tensor. In two limits these give these models reduce to the f(R) models, 
in the metric and in the Palatini formulation. In addition, one may consider coupling 
the scalar to the Gauss-Bonnet invariant. Finally, by considering nonlocal modifications 



By this we mean the manifest form of the field equations. Of course we can reconstruct models 
with identical background expansion to the one implied by the phenomenological Friedmann equation 
involving the integrals (39). 




(39) 




(40) 
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one may avoid the appearance of a finite number of poles into the propagator, whilst 
(multi) scalar- tensor modeles may reflect some properties of these theories. 

The Gauss-Bonnet model can feature a transition from scaling era to dark energy dom- 
ination, but there are problems with consistent solar system tests and new instabilities at 
cosmological scales. The Palatini models may viably reproduce the observations at both 
the cosmological and smaller scales, if one assumes some non-cold properties for the 
dark matter component. The definite predictions of the nonlocal models are presently 
less robustly understood, but further study of such models is worthwhile as they seem 
promising in alleviating the cosmological fine-tuning problems. 
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